Određivanje uopštene kompleksne funkcije
assume(0 < R & 0 < C & 0 < L )
zamena = L == C*R^2
zamena = 
vrednosti = [C == 1, R == 1]
vrednosti = 
Z1 = R + 1/(s*C)
Z1 =

Z2 = 1/(1/(s*L) + 1/R)
Z2 =

H(s) = subs(Z2/(Z1 + Z2), lhs(zamena), rhs(zamena))
H(s) =

simplify(H(s))
ans =

Amplitudska i fazna karakteristika
Hjw(w) = subs(simplify(H(s)), s, 1i*w)
Hjw(w) =

Aw(w) = simplify(abs(Hjw(w)))
Aw(w) =

phi(w) = expand(angle(Hjw(w)))
phi(w) =

Crtanje amplitudske i fazne karakteristike
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
fplot(w, subs(phi(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
Određivanje propusnog opsega
Aref = limit(Aw(w), w, Inf)
Aref =

w3dB = solve(Aw(w) == Aref/sqrt(2), w)
w3dB =

BandPass3dB = [ w3dB, Inf]
BandPass3dB =

wg1 = subs(w3dB(1), lhs(vrednosti), rhs(vrednosti))
wg1 =

Označavanje propusnog opsega u grafiku
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
rectangle('Position',[double(wg1), 0, 4, double(Aref/sqrt(2)) ], ...
'FaceColor',[1 0.9 0.9],...